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Abstract. Recent sampling theorems allow for the recovery of operators with bandlimited 
Kohn-Nirenberg symbols from their response to a single discretely supported identifier signal. 
The available results are inherently non-local. For example, we show that in order to recover a 
bandlimited operator precisely, the identifier cannot decay in time nor in frequency. Moreover, a 
concept of local and discrete representation is missing from the theory. In this paper, we develop 
tools that address these shortcomings. 

We show that to obtain a local approximation of an operator, it is sufficient to test the 
operator on a truncated and mollified delta train, that is, on a compactly supported Schwarz 
class function. To compute the operator numerically, discrete measurements can be obtained 
from the response function which are localized in the sense that a local selection of the values 
yields a local approximation of the operator. Further, we exhibit that the derived mesurements 
allow for signal processing methods based on local features, such as coarse quantization. 

Central to our analysis is to conceptualize the meaning of localization for operators with 
bandlimited Kohn-Nirenberg symbol. 
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1. Introduction 

In communications engineering, the effect of a slowly time-varying communication channel is 
commonly modeled as superposition of translations (time shifts due to multipath propagation) and 
modulations (frequency shifts caused by Doppler effects). In order to recover transmitted signals 
from their channel outputs, precise knowledge of the nature of the channel is required. A common 
procedure for channel identification in this sense is to periodically send short duration test signals. 
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The resulting outputs are then used to estimate channel parameters which allow for an inversion 
of the operator HI H E3 022 [Q [T3] . 

Kailath [14] and Bello [2] analyzed the identifiability of such channels. In mathematical terms, 
these channels are characterized by bandlimited Kohn-Nirenberg symbols and the channel identi- 
fication problem becomes an operator identification problem: can an operator with bandlimited 
Kohn-Nirenberg symbol be identified from the output corresponding to a given test input signal? 

Kozck and Pfander [T5] , and Pfander and Walnut gave mathematical proof of the assertions 
by Kailath and Bello that a suitable test signal exists as long as the band support of the symbol 
of the operator has outer Jordan content less than one. The suggested test signals are periodically 
weighted regularly spaced Dirac-delta distributions. In 22 , Pfander coined the term operator 
sampling as the resulting theory has many direct parallels to the sampling theory for bandlimited 
functions. For example, an operator sampling reconstruction formula was established which gener- 
alizes the reconstruction formula in the classical sampling theorem for bandlimited functions (see 
[22] and Theorem [231 below) . 

The operator sampling results in [15, 25, 22, 24] rely on using test functions as those described 
above. These decay neither in time nor in frequency and cannot be realized in practice. In this 
paper, we show that indeed, for stable identification of operator classes defined by a bandlimitation 
of the Kohn-Nirenberg symbol, test signals that lack decay in time and frequency are necessary. 
When seeking to recover only the operator's action on a time-frequency localized subspace, then 
this ideal but impractical signal can be replaced with a mollified and truncated copy; the test 
signal can thereby be chosen to be a compactly supported Schwartz function as shown below. 

Furthermore, an important difference to the sampling theory for bandlimited functions is that 
the response to a test signal in operator sampling is a square-integrable function rather than a 
discrete set of sample values. Of course, one can obtain a discrete representation using any basis 
or frame of the range space of square integrable functions, but the question remains which of the 
multitude of commonly considered representations allow to recover the operator most efficiently. 
In the case of a bandlimited function, one feature that distinguishes the representation by samples 
is locality: a sample is the function value at a given location; due to the smoothness of bandlimited 
functions it represents the function in the neighborhood of the sampling point. This feature has 
many fundamental approximation theoretic consequences. First, it allows to approximate the 
function in a given region using only samples taken in a fixed-size neighborhood of it. Second, 
locality is a key ingredient for many signal processing methods. Examples are coarse quantization 
schemes such as Sigma-Delta modulation, where the quantization accuracy depends on the good 
control over finite differences between neighboring samples. 

In this paper we develop discrete representations of operators with bandlimited Kohn-Nirenberg 
symbols that, on the one hand, can be computed in a direct and simple way from the output 
corresponding to a test signal and, on the other hand, have locality properties analogous to those 
we appreciate in the classical sampling theory. We work with the same concept of locality as in 
the localized sampling results mentioned above, namely, locality will be defined through the action 
of the operator on time-frequency localized functions. 

Combining the results of our analysis, we show that a suitable local discrete representation 
arises from time-frequency measurements of the output corresponding to a truncated and molli- 
fied weighted sum of Dirac delta distributions. Moreover, these discrete measurements allow for 
arbitrary oversampling rates, which allows to design coarse quantization schemes based on the 
resulting representations. 

The paper is organized as follows. We summarize our main results and put them in the context 
of previous work in Section [2] In Section [3] we recall operator sampling terminology in some 
detail. Section [4] provides results on local approximations of operators; in Section [5] we discuss 
identification using smooth and finite duration test signals. Section [6] uses Gabor frames to derive 
our novel discretization scheme for operators with bandlimited Kohn-Nirenberg symbols, and in 
Section [7] we use the resulting representations to devise a coarse quantization scheme for such 
operators. 
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2. Main Results 



Every bounded linear operator on L? (R) has a formal Kohn-Nirenberg symbol representation 

(2.1) Hf(x)= J o- H (x,0f(0e 2 ™* <% 

where an is a tempered distribution on R 2 , the integral is understood to converge weakly, and 
here and in the following J denotes integration over R. The Fourier transform / is normalized by 

•77(£) = /(0 = / f(t)e- 2 ^dt 

for integrable /. 

The space of bounded operators whose Kohn-Nirenberg symbols are bandlimited to a given 
set M — we will also use the shorthand terminology bandlimited operators — is called operator 
Paley- Wiener spacQ, it is denoted by 

OPW{M) = {H e C(L 2 (R)) : supp T s a a C M}, 

where the symplectic Fourier transform T a is defined densely by 

(2.2) 7>(i, i) = JJ <r(x, 0e- 27ri( ^-*« dx d£ . 

The Kohn-Nirenberg symbol of an L 2 -bounded operator with supp.F s (7i7 compact is bounded. 
In fact, for some A, B > we have, 

(2-3) -4||0ff|U°°(K) < \\H\\c(l 2 (wl)) < B\\a H \\l°°(r) , 

for all H <E OPW(M), where ||-ff ||,c(l 2 (r)) is the operator norm of H (Proposition 14. 1 1 below) . 

Certainly, if we have direct access to an, then some of our approximation theoretic goals can 
be accomplished using classical two-dimensional sampling results applied to an- In the model 
considered here, however, we do not have access to any of the values of the symbol o~h of the 
operator H directly, but we must rely on the operator output Hw resulting from applying H to 
a single test input w. Due to stability consideration, we say that the linear space OPW(M) is 
identifiable by w if for A, B > we have 

(2-4) A\\H\\ c(L 2 m < \\Hw\\ L 2 (m < B\\H\\ ciL 2 {R)) , 

for all H e OPW(M) [15]. "Sampling" the operator means that the identifier w in (|2.4j) is a 
weighted sequence of Dirac delta distributions, that is, 

W = C k S k T, 

where Cfc is an appropriately chosen periodic sequence [THJ [251 122] ■ 

A guiding paradigm in the sampling theory of operators is the direct analogy to sampling of 
bandlimited functions. To illustrate this analogy, we compare the classical sampling theorem (often 
credited to Cauchy, Kotelnikov, Shannon, and Whittaker, among others), Theorem 12. 1[ with the 
corresponding result for operators, Theorem 12.21 [22]. Note that Theorem 12.11 formally follows 
from Theorem 12.21 by choosing the operator H in Theorem 12.21 to be the pointwise multiplication 
operator / i-> a ■ f [2"2"] . 

The engineering intuition underlying sampling theorems is that reducing a function to periodic 
samples at a rate of 1/T samples per unit interval corresponds to a periodization with shift 1/T in 



In general terms, operator Paley- Wiener spaces are denned by requiring its members to have bandlimited Kohn- 
Nirenberg symbol which are in a prescribed weighted and mixed L p space | 22 | . For example, to restrict the attention 
to bandlimited Hilbert-Schmidt operators, we would consider only operators with square integrable symbols. These 
form a subset of the operators considered in this paper. 
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n si 

via a convolution with a low-pass kernel, that is, a function (f> that satisfies 



frequency space [3D]. Thus, as long as TO < 1, a function bandlimited to [— §, §] can be recovered 



(2.5) 0(0 



im if lei < § , 
o, if iei > a? ■ 



If TO = 1, the only such function is <j>(t) = sinc("7rf/T) = S1D ^J^ . For TO < 1, there are many 
such functions; in particular <f> in the Schwartz class is possible. With this notion, the classical 
sampling theorem reads as follows. 

Theorem 2.1. For g e L 2 (R) with suppJg C [-§, §] and TO < 1, we have 
(2.6) 5 (a;) = ^. 9 (nT)0( a; -nT) 

with uniform convergence and convergence in L 2 (JBL). Here, <f> is any low-pass kernel satisfying 
(1231). 



Recall that every operator H on L 2 (M.) is in one-to-one correspondence with its kernel kh, that 
is, for a unique tempered distribution kh, we have Hf(x) = f Kn(x,y) f(y) dy weakly. In the 
following, xa denotes the characteristic function of a set A. 

Theorem 2.2. /|2"2]/ for H : L 2 (R) — > L 2 (R) with <j h G L 2 (M 2 ), supp 7" s ct h C [0,T]x[-f , §], 
and TO < 1, w /lave 

(2.7) K ff (i + t J i) = X[o,n(t) 51 (H^hT)(t + nT)<j>(x-nT), 

n£Z fcGZ 

wt/j convergence in L 2 (R 2 ) and uniform convergence in x for each t. Again, <f> is any low-pass 
kernel satisfying (I2.5[) . 

We point to an important difference between the applicability of Theorems 12.11 and 12.21 in 
Theorem 12. 11 a bandlimitation to a large set [—§,§] can be resolved by choosing a small T; on the 
other side, Theorem 12.21 is not applicable if the bandlimiting set [0,T]x[— ^, has area greater 
than one. Indeed, in [25l [23] the following is shown. 

Theorem 2.3. OPW(M) is identifiable in the sense of (|2.4I) with appropriate w = 5Z ngZ c n^nT 
if M is compact with measure less than 1. If M is open and has area greater than 1, then exists 
no tempered distribution w identifying OPW(M). 

Hence, it is necessary to restrict ourselves to operator Paley- Wiener spaces defined by compact 
sets M with Lebesgue measure one. For such spaces, one can extend Theorem 12. 21 to the following. 

Theorem 2.4. /|24] / Let M be compact with Lebesgue measure less than one. Then exists T, O > 
with TO = j-, L prime, 6 > 0, and an L-periodic sequence (c n ) so that with w = ^2 n c n 8 n T, we 
have for H € OPW{M) 

L-l 

(2.8) K H (x + t,x) = LTY J K* - k 3 T ) ( bji 11 ™^ ~ (*j ~ l) T )<t>(x + (kj - q)T)) e ^ in ' ax , 

j=0 9 GZ 

where r, <p are Schwartz class functions that satisfy 

r(t)£( 7 ) = if (t, 7 ) i (-6, T + S)x (-6, + 5), 

and 

(2.9) £r(t-*T) = l = J2$h -«n). 

fceZ nGZ 

Moreover, (|2.8[) converges in L 2 (R) and uniformly in t. 
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In both sampling scenarios, working with Schwartz class kernels r,<j) is of advantage. Indeed, 
in the classical sampling theorem, the slow decay of the sine function sin x fx in (|2.6p implies 
that a small perturbation of just a few coefficients g(nT) can lead to significant deviations of all 
values g(t) outside of the sampling grid TZ; this includes values achieved at locations far from the 
sampling points nT. Hence to approximately recover the function values locally, that is, on an 
interval [a, 6], it does not suffice to know the function samples in a constant size neighborhood of 
that interval. When working with Schwartz class kernels, in contrast, such a local approximate 
reconstruction is possible; one has 

(2-10) g(x)a £ g(nT) <j>(x - nT), 

n£ [a — d,b+d] 

where the neighborhood size d does not depend on the interval [a, b]. 

A corresponding possibility of using local information for local reconstruction is not given in 
Theorem l2.2l Moreover, the identifier w = X^nez &nT ne hher decays in time or in frequency, clearly 
showing that in practice, this input signal is not usable. However, in the framework of Theorem l2.21 
this is unavoidable, as the following theorem shows. 

Theorem 2.5. If the tempered distribution w identifies OPW{[0, T] x[-fi/2, 0/2]), Tfl > 0, then 
w decays weakly neither in time nor in frequency, that is, we have neither 

(w, ip(- — x)) x ^f° o nor (w, ip(- — £)) o 

for all Schwartz class functions if. 

We address this problem by developing a concept of "local recovery" of an operator, in analogy 
to the local recovery of a function in (I2.10[) . Indeed, the key to most results presented in this paper 
is to aim only for the recovery of the operator restricted to a set of functions "localized" on a 
prescribed set S in the time-frequency plane. This is indeed reasonable in communications where 
band and time constraints on transmitted signals are frequently present. In [13 , for example, 
operators that map bandlimited input signals to finite duration output signals are considered. 
Bivariate Fourier series expansions of such an operator's compactly supported Kohn-Nirenberg 
symbol allows the authors to discretize the a-priori continuous input-output relations (|2.ip and 

The definition of functions localized in time and frequency is based on Gabor frames. Their 
definition involves translation and modulation operators, 

71/ :/-►/(• - t) and M u : f ^ e 2 ^f. 

These operators are unitary on L 2 (M) and isomorphisms on all function and distribution spaces 
considered in this paper. 

For any g £ i 2 (M) and a, b > 0, we say that the Gabor system 

(g,aZ x bZ) = {Tk a Me b g} k ,eez 
is a tight frame for L 2 (M. d ) if for some A > 0, we have 

f = Aj2 (f,T ka M ib g) TkaMetg 

for all / £ L 2 (R d ). Each coefficient in this expansion can be interpreted to reflect the local behavior 
of the function near the indexing point in time-frequency space. Hence, a natural way to define 
time-frequency localized functions is that all but certain expansion coefficients are small. 

Definition 2.6. Let (g, aZ x bZ), g £ <S(R), be a tight frame for L 2 (M.) with frame bound 1. We 
say that f £ L 2 (R) is e-time-frequency localized on the set S if 

E \UMibT ka g)\ 2 >(l-e 2 ) J2 \(f,M eb T ka g)\ 2 . 

(ka,£b)eS (ka,eb)£R 2 
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Our first result states that a sufficient condition for two operators to approximately agree on 
functions e-time-frequency localized on a set S is that their Kohn-Nirenberg symbols almost agree 
on a neighborhood of S. Below, B(r) denotes the Euclidean unit ball with radius r and center 0; 
the dimension will always be clear from the context. For brevity of notation, we set S — B(r) — 
(S c + B(r)) c for S C R 2 . 

Theorem 2.7. Fix M compact with fi(M) < 1 and let (g, aZ x bZ), g £ <S(R) ; be a tight frame for 
L 2 (R) with frame bound 1. Then exists G > and a strictly monotone function d : (0, 1) — > M + , 
lim^o d(e) = +oo, with the property that if H.H £ OPW(M) satisfy on a set SCl 2 the bounds 

\\o-h\\l~(rz), ||o#||z,oo (R 2) < [i and \\a H - 0fi\\L°°(S) < ^Mi 

then 

\\Hf-Hf\\„ (R) <Cen\\f\\ L , {R) 

for all f £ L 2 (R) that are e -time- frequency localized on S — B(d(e)) in the sense of Definition \2.b\ 

This observation is a key ingredient in the proof of our next main result. It concerns truncated 
and mollified versions of the identifier J2 n c n5 n T and provides localized versions of Theorems 12.21 
and|U For S = K 2 , it reduces to Theorems I2T21 and [2^1 

Theorem 2.8. Fix M compact with fi(M) < 1 and let (g, aZ x b%), g £ <S(R) ; be a tight frame for 
L 2 (M) with frame bound 1. Choose S > such that n(M + [—3(5, 35] 2 ) < 1. Then exists C > and 
a strictly monotone function d : (0, 1) — > R + , lim e _>.o d(e) = +oo, with the following property: 

Let H £ OPW(M) satisfy the bound W&hWl^Cr 2 ) ^ M an d let S £ Ii x I 2 £ M 2 , where 1\ and 
I2 may coincide with R. Furthermore, let 

where the tempered distribution if is chosen such that ip > and <p = 1 on I2 and define for 
H £ OPW(M) the operator H via 

L-l 

(2.11) n R (x + t,x) = LT^2 r(t - k 3 T) ( b jq Hw(t - (kj - q)T)(f)(x + - q)T)) e 2 ^ na: , 

i=o gez 

where r,(f> are Schwartz class functions defined as in Theorem \2.4\ but for the above choice of 5. 
For rectangular bandlimitation domains M = [0, T] x [-S W] one can choose the identifier 
^2 tp(- — nT) and define H via the formula 

nTe/i 

Kjj(i + (,i)=T^(ff f{- - nT)){t + nT) 4>{x - nT) . 
In both cases one has 

\\Hf-Hf\\ LHR) <Ce^\\f\\ L 2 {m 

for all f £ i 2 (R) that are e -time- frequency localized on S — B(d(e)) in the sense of Definition ] 2. b\ 

Note that this theorem is completely analogous to the condition (|2.10[) for localized function 
sampling. Due to the two-dimensional nature of the operator, however, localization is an issue 
in both time (restricting to a finite number of deltas) and frequency (replacing the deltas by 
approximate identities). If one is interested in localization only in time or only in frequency, one 
can choose one of the Ii to be R and thus consider 

w = 2_j c ^nT or w = } j c n ip{- - nT), 

nT£h n 

and with (c„) = 1 in case of rectangular domains M. 
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An additional important structural difference between classical sampling and operator sampling 
remains: in Theorems 12.21 and 12. 81 the reconstruction formulas (|2.7p and (|2.8I) involve as "coeffi- 
cients" functions, not scalars. Among the many possibilities to discretely represent the operator's 
response to the identifier w, we consider Gabor representations of this sample function. A time- 
frequency localized subset of the coefficients will then yield a corresponding local approximation of 
the operator. Theorem 12.91 below establishes a reconstruction formula based on Gabor coefficients 
that allows for the exact recovery of the operator; Theorem 12.101 shows that a local subset of the 
coefficients yields a local approximation of the operator. Again, one obtains considerably simpler 
formulas for rectangular domains, but for reasons of brevity, we focus on the comprehensive setup 
of arbitrary domains. 

For a Schwartz class function 4> and a tempered distribution / on K we call 

(2.12) Vtf(x,t) = (f,M t %<f>), i^ei, 

the short-time Fourier transform of / with respect to the window function (f>. Throughout this 
paper, all pairings (•, •) are taken to be linear in the first component and antilinear in the second. 

Theorem 2.9. For M compact with fi(M) < 1 exists L prime, 8 > 0, T, Q > with Tfl = l/L, 
an L-periodic sequence (c n ), and a sequence {-Bj g }j=o....,L-i.ij(EZ (which is L periodic in q and 
depends only on the sequence (c n ) ), such that with w = ^2 n c n S n T we have for H € OPW(M) 

(2.13) 

where 

*m,e = Y, B ™ 0(HZ- fcj-mL/A)r) (Hw, T qT M mL/ p 2 r), 

and r, <j> are Schwartz class functions such that r and <f> are real valued and satisf^ 

(2.14) r(t) = if t<£ {-8, S + T), £( 7 ) = if 7 $ (-5 - 0/2, 8 + 0/2), 
and 

(2.15) £ \r(t + kT)\ 2 = 1 = E l?(7 + ^)| 2 , 

fcSZ nGZ 

with oversampling rates fa > 1 + 28 /T and fa > 1 + 2<5/O0 

Observe that the reconstruction formulas given in Theorems 12.41 and 12.81 require r and <f> to 
generate partitions of unity (|2.15[) . while (|2 .9[) above requires that their modulus squared form 
partitions of unity. 

Theorem 2.10. Fix M compact with fJ-(M) < 1, let T, O, L and w, r, <f> be defined in Theorem \2.4[ 

and let (g,aZ x UL), g G 5(R), be a tight frame for L 2 (M) with frame bound 1. Then exists C > 
and a strictly monotone function d : (0, 1) — > R + , lim e _j.o d(e) = +00, with the following property: 
Let H e OPW(M) satisfy the bound ||<xh-||i,°°(k 2 ) < I 1 , fix a not necessarily bounded set S C 
I\ x I2 in R 2 , choose ip and w as in Theorem \2.8\ and define the operator H via its symbol 

PlP2 3=0 ( m LT/p u lLSl/p a )eS Pl ?2 



2 For example, we can choose r = X[o,T) *fSt where <ps is an approximate identity, that is, a non-negative function 
with ipg G <S(R), snppips Q [ — <5/2, <5/2] , and f <ps = 1. 

3 Then the Gabor systems {r fej ; = T tT M(/feT rltia, {T„nM ro /ftn?}m,neZ, and = 
T( mT L/l3- L ,tLn)M(nn,/j32,kT}}m,n,k,t& are tight Gabor frames with A = fa/T, A = fii/Q, and A = /3i/3 2 /(TC) = 
plfoL, respectively, whenever /3 2 > 1 + 28/T and fix > 1 + 28/Q. 
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where 

° ( £i = E B ™ - k i ~ mL /Pi) T ) (Hw, T&MtnL/h r). 

gGZ 

Then H satisfies 

\\Hf -H.f\\ LHm <Ce/i|l/IU»(R) 

/or / € L 2 (R) which are e-time-frequency localized on S — B(D(e)) with respect to (g,aZ x 6Z) 
m i/ie sense of Definition \2. 6\ 

The discrete representations introduced in Theorems I2.9l and l2.10l resolve a fundamental concep- 
tual difference between classical sampling and operator sampling. In contrast to classical sampling, 
which yields a set of separate function values, the contributions of the different Dirac-deltas in the 
operator sampling formula are combined in a single function and cannot easily be separated. Hence, 
while choosing a higher sampling rate in the function case yields more information, in the oper- 
ator case, this additional information is mixed in an inseparable way. These aliasing effects |15) 
make it impossible to obtain redundant representations merely by oversampling in Theorem 1 2 . 2 1 or 
Theorem 12.41 In reconstruction formula (|2.13[) . however, the oversampling parameters can be 
chosen arbitrarily, allowing for representations of arbitrarily large redundancy. 

The interplay of large redundancy and good local representation properties of the discrete co- 
efficients allows for coarse quantization methods to be applied to operators. In the mathematical 
literature, the most common scenarios for such methods deal with frames in 1" [3J HI [T7] or the 
space of bounded bandlimited functions on R [5J [T21 [7] . The underlying idea is to sample at a 
high rate and use the resulting redundancy to lower the number of bits needed to represent each 
sample. The following establishes a corresponding result for bandlimited operators. 

Theorem 2.11. Fix a compact set M with fJ-(M) < 1 and the associated parameters L, Q, and 
T as in Theorem \2.9l Then there are constants c, C such that for each H G OPW(M) satisfying 
\\o~h ||l°°(k 2 ) < c and all oversampling rates Pi, 02 > 1, one can compute quantized values q^' \ G 
{±1, ±3} directly from Hw with w = CkSkT via a recursive procedure. These quantized values 
give rise to an approximate reconstruction of H through the operator H with symbol 

a(x,0 = ^j:e- 2 ^^^e^^ £ V,r(x - + kj )T, £-(f +nj p), 

which satisfies 

\\Hf-Hf\\ L2{m <JL||/|| La(R) 
P\Pi 

for all f G L 2 (R). 

Combining all these results, we obtain a result about local reconstruction from local quantized 
values resulting from a localized identifier. 

Corollary 2.12. Fix a compact set M with fJ.(M) < 1 and the associated parameters L, Vt, and T 
as in Theorem \2.9\ Then there are constants c, C , and a strictly monotone function d : (0, 1) — > 
R + , lirn e -H> d(e) = +oo such that for each H € OPW(M) satisfying \\cth\\l°°(m 2 ) < c an d °^ 
oversampling rates /3i,/?2 > 1? one can compute quantized values q^i € {±1,±3} from Hw, with 
w as in Theorem \2.8\ via a recursive procedure. Then using the finitely many quantized values with 
indices (m,£) satisfying (mLT '/ ' (3\, £LQ / ' P2) G S, one obtains a local approximate reconstruction of 
H on S — d(e) through the operator H with symbol 

PlP2 j=0 (mLT/t3 1 ,lLn/f3 2 )eS Pl P2 
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This operator satisfies 



~ C 
\Hf - Hf\\ L 2 W <—— ||/|| i2(R) +Cen\\f\\ L 2 W 

P1P2 



for all f € L 2 (R) which are e-time- frequency localized on S — B(D(e)) with respect to (g,oJL x 6Z) 
in the sense of DefjnitionW.~ 



3. Bandlimited pseudodifferential operators and operator Paley- Wiener spaces 
It is well known that every bounded linear operator H : SiW 1 ) — > S^R^) is of the form 

Hf(x) = J K (x,t)f(t)dt 

for some k € <S'(R 2d ), where 6>(R d ) is the Schwartz space, and S'(M. d ) is its dual, the space of 
tempered distributions |10j . This integral representation is understood in the weak sense, that is, 

(Hf,g) = (kJ® 9 ) 

for all /, g G S(R d ), where f®g{x, y) — f(x)g(y) and ( ■ , • ) is the sesquilinear pairing between S 
and S' functions. Each such operator has a spreading function representation 



(3.1) Hf = J J v(t,l)M y T t f dtdrt, 

a time-varying impulse response representation 

Hf(x) = J h(x,t)f(x-t)dt, 

and a Kohn-Nirenberg symbol representation 

Hf(x) = J a{x,i)m^ ixi 

We write H a and an, gn, hh when it is necessary to emphasize the correspondence between H 
and a, 77, k. The symbols a and 77 are related via the symplectic Fourier transform J- s defined in 
(|2.2p . that is, a = T s g- For convenience, we use the symbol rj H (t,-f) = e 27Tnt gii{t, 7), and denote 
its symplectic Fourier transform by 



The relationship between an and er h is given by 
(3.2) oh* = cth, 

where -ff * denotes the adjoint of H . Indeed, we have for Schwartz functions /, g 
(H*f,g) = (f,Hg) = (f(x), ( [ r) B (t,v)M u T t g(x) dt dv) 



m(t,v){f,M„T t g} dtdv 



VH(t,v)(T-tM-„f,g) dtdv 



VH(-t, -v)TtMvf dtdv, g) 



n H (-t,~v)e- 2 ^ MJTtf 'dtdv, g) 



and conclude using a density argument that t)h* (t, v) = T)h{~ t, —v) e 2mtv . Hence, 
oh. (x, = F.w {*,( ) = jj nH(-t,-w)e- 2 '"- €.-*"•• e 2 "" dt dv 

= J J r] H (t, v)e 2 ^ itv e ~ 2nit S e 2 * ixu dt dv 

To prove our results, we shall frequently transition from a to er. This does not cause a problem 
in our analysis since (|2.3|) combined with \\H ||z:(i 2 (R)) = \\H* \\c(l 2 (Wl)) shows that for M C R 2 
compact exist A, B > with 

(3-3) ^lklli«>(H 2 ) < lklk~(R 2 ) < B\\cr\\ L o 0{R 2 ) 

for all H a e OPW(M). 

4. Local approximation of bandlimited operators 

In this section we show that a local approximation of an operator's symbol always yield a local 
approximation of the operator in the sense of Definition ^. 61 The given results are of general interest 
and will be stated in more general terms than other results in this paper. This does not increase the 
difficulty of proof, but necessitates to recall additional terminology from time-frequency analysis. 

For that, recall that for any full rank lattice A = AZ 2d C M 2<i , det A ^ 0, £ P (A) denotes the set 
of sequences (ca)asA for which 

iicii^(A) = (£M p ) 1/P<o °- 

AeA 

A time-frequency shift by A = (t, v) G A is denoted by 7r(A) = A4„Tt and in the following we will 
consider Gabor systems of the form (g, A) = {7r(A)g}AeA- 

Among the many equivalent definitions of modulation spaces, we choose the following. Let 
go(x) = e" x " and 1 < p < oo. Then 

(4.1) M p (R d ) = {/ G S'(R d ) : ||/|U (Rd) = \\{{fM^9o))xh Hi z^) < ™} 

(see, for example, [TTJ [H]). In the following we shall use the fact that whenever (g, A) is a tight 
i 2 -Gabor frame (see below for a precise definition) with g £ Af 1 (M <1 ) then replacing the L 2 -Gabor 
frame (go, \l? d ) in (|4.ip with (g, A) leads to an equivalent norm on M p (R d ) [11]. That is, there 
exist positive constants A and B with 

(4-2) Mf\\ P u, m < E K/^( A ).9>I P < S II/IIm,(^)' f £ MP ( Rd ) 

AeA 

if 1 < p < oo and 

A\\f\\Moo {m <su V \(f,n(X)g)\<B\\f\\ Mxm , f G M°°(W. d ) 
AeA 

if p — oo. In either case, we call (g, A) an £ p -frame with lower frame bound A and upper frame 
bound B. If we can choose A = B in case of p = 2 then we call (g, A) a tight Gabor frame. 
The norm equivalence (|2.3p follows from the following result since M 2 (R) = L 2 (M.). 

Theorem 4.1. Let 1 < p < oo and M compact. Then exist positiv constants A = A(M,p) and 
B = B(M,p) with 

A \\<th\\l-(F) ^ W H Wc(M P (R)) < B \\a H \\ L oo m , H G OPW(M). 
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Proof. Theorem 2.7 in [35] (see for example the proof of Theorem 3.3 in [32]) provides C = C(M,p) 
with 

||#/I|m»(R) < C \\o-h\\l°°(r 2 ) ||/||mp(r) 

for all H e OPW(M). This establishes the existence of B = B(M,p) above. 

In addition, we shall use the following facts. In [9] [IT] it is shown that the operator norm 
of an operator mapping the modulation space M 1 (M) into its dual M°°(R) is equivalent to the 
M°°(R 2 ) norm of its kernel k, which can easily shown to be equivalent to the M°°(R 2 ) norm of 
the time- varying impulse response h. Moreover, we use the fact that M°°(R 2 ) is invariant under 
Fourier transforms (in some or all variables) and that the M°°(R 2 ) norm can be replaced by the 
I/°°(R 2 ) norm if we restrict ourselves to functions bandlimited to a fixed set M [TH1[22]. Last but 
not least, we use that the identity map embedding M P (R) into M 9 (R), p < q, is bounded. 

Writing < to express that A < CB for some constant C depending only on the support M and 
Ax5to denote equivalence in norms, i.e., A < B and B < A, we obtain for all H £ OPW(M) 

\\&h\\l°°(R 2 ) ~ llcffll jW"e>°(R 2 ) X ll^-H ||m°°(B 2 ) X \\kh\\m-=°(R 2 ) 
x ||-^IU(M 1 (M),M~(R)) i$ \\H\\c(Mp(W)) 

and the result follows. □ 



We proceed to prove the following generalization of Theorem 12.71 Indeed, the earlier stated 
result follows again from the fact that L 2 (R) = Af 2 (R). We focus on the case of arbitrary domains; 
a simpler proof for rectangular domains can be obtained using Theorem 12 . 2 1 instead of Theorem l2.4l 

Theorem 4.2. Fix M compact and p € [1, oo] . Let (g, A), g £ M 1 (R), be a tight frame for L 2 (R) 
with frame constant 1. Then exists a constant C and a strictly monotone function d : (0, 1) — > R + , 
lim £ _>.o d{e) = +oo, with the property that if H £ OPW(M) satisfies 



the 



\o-h\\l^(m 2 ) < M and \Wh\\l°°(s) < 



\\Hf\\Mp(R) < C e/j.\\f\\ MP / 



for all f £ M P (R) time-frequency localized on S — £>(d(e)) = yS c + B{d{e))j in the sense that, 
for p < oo, 

E K/,7r(A) fl >| p >(l-e*) EK/'^Wr 

AeAn(S-B(d(e))) AeA 

or, for p = go, 

sup{K/,7r(A)ff)|, AGAn(5-B(d(e)))}>(l-e) B up{K/,7r(A)ff)|, AeA}. 

Proof. Step 1. Preliminary observations and choice of auxiliary objects. Choose a nonnegative 
<f> € S(R 2 ) with / 4>(x) dx = l and supp (j> C [-|, ±] 2 . Recall that 

A 1 = {fi e I 2 : e 2^<M:A> = j for aU A G A | 

is called dual lattice of the lattice A in R 2 . Let A be a lattice containing A with the property that 
there exists a compact and convex fundamental domain D of A -1 which contains M + [— i, i] 2 . Set 

o-p = Wxd * '/'II^V 2 ) ^(xd * 4>) 
and, using the sampling theorem for lattices in R n [2TJQJ], we obtain for all H £ OPW(M) 

o-h = 2J ct h(A) 7a o-p 
AeA 

and hence 

(4.3) H = E^h(A)^(A)Ptt(A)*. 

AeA 

11 



As explained above, the fact that (<?,A) is a Gabor frame in L 2 (R) with g E M 1 (M.), implies 
that it is also an £ p -frame for Af p (R) and there exists C\, C2 > with 

(4.4) \\f\\ M p(R) < Ci ||{(/,7r(%)} AeA ||^ (A) < C\C 2 \\f\\ Mm , f G AP(R). 

As the synthesis map is the adjoint of the analysis map, we also have 

(4-5) || c A7r(*)ff|| M »(K) - ° 2 II{ c a}aga|Up(a). 

AeA 

For some /ii, /12, • ■ • , Mm we have A = U"=i(^ + w)i where n depends only on Af and (3, A). It is 
easily seen that (g, A+/i^), ^ = 1, . . . , n, also satisfies (|4.4[) and (14.51) . Setting 3 = rT 1 ! 2 g e M 1 (M), 
we conclude that the Gabor system (g, A) is a tight frame for L 2 (R) with frame bounds equal 1 
and an £ p -frame with for AP(R) with 

II/IUp(r) < Cin*~» ||{</,7r(A)5)} Xe xll ff( A) 

< C 2 np"i ||/|U, (R ) = CiC 2 ||/|| MP ( R ), / G AP(R). 

We claim that 

{{P*(\)g,w(\)g)} GZ^AxA). 

To see this, recall that er P G S(R 2 ) C Af X (R 2 ), and, hence, a P given by a P (x, £) e 2mx ^ is in M 1 (R 2 ) 
as e 2lx%xS - is a Fourier multiplier and hence also a time multiplier for M 1 (R 2 ) (Theorem 11, [5]). A 
direct computation implies that for A = (t, v) and A = (i, v) we have 



\{P7r(t,u)g,7r(t,V)g)\ = a P (x,0 e 2 ^ M v T t g{§Mi;T~g{x) d^dx 



a P (x, e 2 "^ M- t %g(OM v Ttf(x) d^dx 
= \ (opM(vA)T { t tV) g®f)\. 

Equation f|4.2[) implies that the right hand side is in I 1 (A x A) since a~p € Af 1 (R 2 ) and (<?<8>g, A x A) 

is a Gabor frame with window g<2)g in M 1 (R 2 ). 
Fix e > and choose d(e) > so that 



J2 E (^(A)g,7r(A)g) 

AeA AeAnB(d(e)) c 



< e. 



Now, set A(X, A) = (Pn(\)g, Tr(X)g) if A € A n B(d(e)) c and else. 

Step 2. Decomposing Hf as Hf = i7 in / in + H out f in + Hf out . 
Set 

A io = An(S-B(rf(e))), A to = An(S-B(d(f))), A out = A\A in , A out = A\A in . 

Let 

fin = E ^' 7r ( A )ff) 7F ( A )-9 = E CA 7r ( A )3' /° ut = / ~ 

A e A in AeA in 
where c\ = yjn (/, 7r(A)g) if A <E A and else. Similarly, inspired by (|4.3|) . we set for H e OPW(M) 

H in = £ ^ff(A)7r(A)P7r(A)*, H out = H - H m 
xeAns 



and note that H m ,H out G OPVF(L> + ±] 2 ). 
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Step 3. Bounding ||-ff ut/in||Mp(R) ■ We use the separation of A; n and AP\S C by d(e) to compute 
I (^out/in, tt(A)5)| = |< E <T H (v)-K(v)P-K{vy E C\ir(\)g, 7r(X)g) 

< E \oh(v)\ E |cA||<7r(i/)P7r(i/)*7r(A)flf,7r(A)flf> 
!/6AnS c AeA in 

< E I^HI E |cA||(P7r(A-^)5,7r(A-^)5) 
velvnS" AeA in 

< E I^hMI E \cx\A{\-v,\-v) 

^eAnS c AeAi„ 

< iio-ffiu~(K 2 ) E E I ca ' ^( a ^ v > a - 

v£A AeA 

For every sequence {c^} € ^ 9 (A), 1/p + 1/g = 1, we conclude 

| ({ (-^out/in, 7r(A)g)}3; eX , {<%} XgX ) | 

< WvhWl-p.') E E E l c A|^(A-f,A-i/)|d x | 
AeA i/eA AeA 



= iiwiiio.^i) E E E I ca +-i A ( A ' A ) i d A+J 

AeA AeA i/eA 

< Ik// IIl°°(k 2 )I|{ca}|| £ p(X) \\{dx}\\ EE^ A ) 



AeA AeA 

and 

\\H out fm\\Mp(W) < ri*~* Ci\\{\(H oui f in ,w(X)g)\}\\ ep( j^ 

< n^p C^lcrHW^^WlcxjW^^ E E A ( A ' A ) 

AeA AeA 

<n'~ Ci/z||{n5 (/,7r(A) 5 )}|| ^p(A) e 

< n 1_ p Ci C 2 e/j, ll/H Mp(a)- 

Step ^. Bounding \\H f out \\ MP ( R) . 



By Proposition 14. II we have 

||-ff/out|Ufp(B) < B(M,p) ||cr//||L~(R2)||/out||Mp(B)- 

By hypothesis, for p < oo we have 

ll/out|lMP(E) = II E (/' 7r ( A )3) 7r ( A )5llMp(R) 

AeA out 

<cf E Ki>(%>l p 

AeA ou , 

^cfe^EK/^wr 

AeA 
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and for p — oo we have 

||/out||M°°(R) = II X! (/' 7r (^)ff) 7r (-^)5lU/-(R) 
AeAout 

<C 2 |K/ l7 r(A) 5 )|Uoo (Aout) 
<C 2 e||(/,7r(%)|U» (A) 

< cf e II/IIm°°(r) ■ 

We conclude 

||#7out||Mp(R) < B(M,p)C 2 e \\<TH\\L°°\\f\\MpQg.) < £(M,p)C 2 e/i||/||M»(R)- 

Step 5. Bounding \ \ H - m fi n ||mp(R)- Since up e 5(]R 2 ), the operator 

£°°(A) L°°(R 2 ), {c A } ^ £ c A T\(7p 

is bounded, say with operator norm bound C3. Then, Proposition 14. 1 1 implies 

1 1 -ffin /in 1 1 Mp (R) < B(D+[~ , i] 2 , p) || <7fl, n || L =o (R 2) || f in || M p(R) 

< B{D+[-\, i] 2 ,p) C 3 || W(A)}|| £=0(AnS) (l + e)||/|| M , (ffi ) 
<2B{D+[-\,\]\p)C i e t i\\f\\ MHm . 

Since all constants are independent of e, /1, H, and /, we summarize 

H-ff/H Afp(R) = ||-ffin/in + -^out/in + -ff /out || Aff(E) < Ce ||/|Ufp(R) • 

□ 

5. Operator identification using localized identifiers 

This section analyzes identifiers that are localized in time and frequency. Theorem 12.51 shows 
that such functions cannot serve as an identifier for the complete operator Paley- Wiener space as 
a whole. 

Proof of Theorem \2.5[ Let r =/= be a Schwartz function with suppr C [0, T] and (j) =/= 
be a Schwartz function with supp0 C [—Q/2,Q/2]. Let H n be defined via its kernel K n (x,y) = 
<j>(x - n)r{x - y), so /i„(ie, t) = 0(x - n)r(t) and r) n (t, u) = J h n (x, t)e~ 2lxixv dx = r(t)e 2 * inv $(y), 
so H n G OiW([0,T]x[-n/2,n/2]) with ||a ff JUoc (B2) = ||f|| L oo (R) ||0|| L <» (R) . 

If tu identifies OPW([0,T]x[-fl/2,n/2]), then by definition H n w G L 2 (M). Then 

J \H n w{x)\ 2 dx = J \(n n (x,y), w{y)) y \ 2 dx 

= J \^(x-n)\ 2 \(r(x-y),w(y)) y \ 2 dx. 

Clearly, (j-(x — y), w(y) ) y x ~^f° Q would imply ||.ffnHI.L 2 (R) "^-T and contradict identifiability 
(EHJ) since by jOJ we have ||i?„|k(L=(E)) > IIl-(r 2 ) = ^||^||i=»(is) ||</>IU°°(r) for all neZ. 

To show that an identifier w cannot decay in frequency, we choose H n € OPW([0, T]x[— ^, ^]) 

to have spreading functions r) n (t,v) = r(t)e 2mnt 4>(iy)e~ 27! " ltL/ . Let j be a Schwartz function and 
compute using Fubini's Theorem and, for notational simplicity, using bilinear pairings in place of 
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sesquilinear ones, 

(H n w^),9(x)) x = (vn(t^),(^ ixv w(x-t),g(x)) x ) tv 

= (r(t)e 2 * int j(v), (e^-^wix ~ t), 9 {x)) x ) tu 
= (r(t)e 2 * mt w(x - t) g(x), e 2 ^^),)^ 
= {{r{t)e 2 ™\w{x-t)4>{x-t)) t ,g{x)) x 
= ((?(£ -n),e- 2 ™tw *m)t,g(x)) x 
= {r(t-n)w*j(0,g(0) f - 

Hence, 

\\H n w\\ 2 L ,, R) = \\G\\ 2 = I \?(t-n)\ 2 \(w(Z-v),$(v)) v \ 2 dt, 



and we can conclude as above. □ 
We proceed by showing that local identification of operators is possible with identifiers localized 
both in time and frequency, Theorem 12.81 

Proof of Theorem \2.8\ The proof proceeds in two steps. First we show that replacing each 
Dirac-delta by a suitable smoothed out version locally introduces only a small error and identifi- 
cation using the resulting smooth identifier can be interpreted as sampling a modified bandlimited 
operator. Second we show that reducing to a finite number of samples also locally yields only 
a small error. Applying this to the modified operator arising in the first part proves that both 
reductions together also yield only a small error. 

For the first part, choose <p € S with supp<y5 C [—5,5], ||^||z,°°(r) = 1, and \@(Q — 1 1 < e for 
£ G I2 • Define C v : f f * ip and set He = H o C v . Observe that 

He fix) = II \ H {t,v)e 2 ™ v f*y{x-t)dtdv 

m {t, v)e 2 ™ u f{x - t - y)<p(y) dy dt dv 
m (t - y,v)e 2 « ixv f{x - t)<p(y) dydtdv 
( / r, H (t-y,vMy)dy) e 2 ^ 'j %x ~ t)dt dv , 

that is, r)H c (ti v ) — Vh(', v) * ip(t) and 

supp?7Hc Q suppi] H + [-5, 5]x{0}. 

We can apply Theorem 12.41 for the operator He with Mi := M + [— 5, 5]x{0} in place of M. As 
by assumption M x + [— 5, 5] 2 still has measure less than one, this can be done with 5, r and <j> as 
given in the theorem. Defining w\ := ip* w, we obtain 

L-l 

(5.1) K Hc (x + t,x)=LT^2 r (* - k 3 T ) ( h iiHw x (t - (kj - q)T)<j){x + (kj - g)T)) e 2 ™^ 

j=0 ?GZ 

Observe that 

<?H G (x,£,) = T a J] Hc (x, £) =<r H (x,0 

and, by hypothesis, we have \\cr Hc \\ L ^(u 2 ) < lk-ff||i~(R 2 ) < /•* and \\a H - o Hc \\l°°(s) < 

Note that for l\ = M, (|5.ip agrees with (|2.11[) and we have He — H, so this establishes the 
result. 

For the second part, let us assume S C I\ x K and Mi C [c, d] x R. Let ip £ S(M) be 
nonnegative and satisfy J2n^( x ~ ~ 1 anc ^ suppV' C [-1/T, 1/T]. Such a function can be 
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obtained by choosing an arbitrary bandlimited, nonnegative ipo € <-> with H^oIIl 1 = 1 and defining 
i> = X[o,T] * i>a- 

Set Pa{x) = J2nTeA ip( x ~ n T), so P[_jv,at] 1 an d P[_Ar.Af]c — > uniformly on compact subsets 
as N — >• oo. Moreover, |Pa(x)| < 1 for all A. Choose N(e) so that | Pf 1 _ ) _[_jv(e),Ar(e)](^) — 1| < e for 
x e i"i + [c, d] and choose R(e) with 



(5- 2 ) 51 \\P[i 1 +[-N(e),N(e)](x)V r r(x-q,Q\\ L i {R 2 ) <e(l-e)D. 

9T^Ji + [-fl(e),fl(e)] 



where the nature of I? is derived by the computations below. The existence of such R(e) follows 
from the fact that J 3 / 1 +[-A'(e) l Ar( e )](a:) and V^-r decay faster than any polynomial. 

Let w 2 = J2kTei 1 +[-R(e) RU)]+[-5 T+8] c k$kT and observe that H as defined in the theorem 
satisfies 



h^(x + t,t) = Kg(x + t, x) 
L-l 

= LT^2r(t - k j T){j^b jq H c w 2 {t - {kj - q)T)^{x + {kj - <z)T)) e 2 ^ Qx . 



Since M\ C [c, d] x R, we have supp Hc5 y C [c + y, d + y], and therefore, 



H c w{x) = Hc^c k S kT (x) = H c ^ c k 5 kT (x) = H c w 2 (x), 

fceZ fere/i + [-fl(e),-R(e)] + [-T-5,«] + [c,d] 

=Ii + [-R(e), R(e)} + [-T -5,5}. 



Note that H e OPW(M 2 ), where M 2 = Ah + [-5,5] 2 (for details, see, for example, [23]). As 
M2 + [— <5, £] 2 still has measure less than one, this implies that we can apply Theorem 12.41 again 
with the same 5. We obtain 



h Ho (x + t,t) — hg(x + t,t) 

L-l 

=LT J2 r(t - kjT) ( b jq H c (w - w 2 ) (t - {k 3 - q)T) <f>{x + {kj - q)T)) e 



3=0 q£Z 
L-l 

=iT^r(t-fcjT)( J2 b jq H c (w-w 2 )(t-(k j - q)T) ^(x + (k 3 - q)T)} e 2m ^ x 

3=0 qT^K-it-kjT) 
L-l 

--LT r(t - kjT) ( Yl b jq H c (w - w 2 ) (t - (kj - q)T) cj>(x + [kj - g)T)) e 2 ™^. 

j=0 gT^/ 1 + [-_R(e),_R(e)] 
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Setting K = K c + [—5,T + S] and using that (a Hc { x ,0 ~ ^(x,^)) Pj 1+ [_jv( e ),iV(e)]( :c ) is ban- 
dlimited to M + {0} x [- 1/T, 1/T] ) , we compute 

\Wh c - 0"jj|U°°(S) < 1/(1 - e) \\((*Hc( x >0 - VH( X '0) P Ii + [-N(e),N(e)](x)h<>°(M?) 

~ 1/(1 - e) IKo-ffc^.O ~ o-#(^)6)- p ri+[-iV(e),Ar(6)](a;)||M~(M2) 

- 1/(1 - e) IIOffc(M) ~ ^( X '*))- P Ji + [-Ar(6),JV(e)](^)||M«'(R=) 

L-l 

x 1/(1 - e ) ||iT P h+[ - N{ e),N(*)]( x ) E r ^ - fc ^ T ) e 2 ""^^-*) 

^ & i9 #c(™ - - (% - g)T) <f>(x-t + (k 3 - 9 )T) 

8T^Xi + [-Ji(e),fl(e)] 
L-l 



M°°(R 2 ) 



< LT/(1 - e) £ E 

j=0 gT^i + l-fite),^)] 



L-l 



< LT/(1 - e) £ E 

j=0 gT^/i + [-/i(e),ii(£)] 



Pz 1+h ;v( e ),Ar(e)](z) r(t - fc,T) e 2-%o(x- t ) 
b jq H c {w - w 2 ){t - (kj - q)T) 4>(x — t+ (kj - q)T) 
H c {w ~ w 2 ){t - (kj - q)T) 



M°°(R 2 ) 



M°°(R 2 ) 



< \\Hc\\c{M^{K)) \\W ~ W 2 \\m 



P h+[ -N(e),N(*)]( x ) r (t ~ k J T ) e 2 ™* n( - x -V b jq <j>{x-t + (kj - q)T) 
LT 



M 1 (R 2 ) 



1-e 



L-l 



M 1 (R 2 ) 



E E ^ P /l + [-W(e),N(e)]( :Z: ) r 'W0(2 ; - i -'? T ) 

where we used the invariance of the M°° and M 1 norm under translation and modulation and, 
for the last inequality, Theorem 14.11 - noting that, for functions constant in one of the coordinate 
directions, the M°° (R) and M°° (R 2 ) norms agree. The second to last inequality is based on M 1 (R 2 ) 
being a Banach algebra, namely on \\gig 2 \\ a/i(k 2 ) < ||.9i||Af 1 (R 2 )ll.92|Uf 1 (K 2 ) for 9i,9i e M 1 (R 2 ). 
Indeed, for / e AP°(R 2 ) and g € Af 1 (R 2 ), we have 

ll/.9lU/-(R 2 ) = _ sup \{fg,f)\= ^ sup \(fjg)\< _ sup H/IImo 

II/IIm1(K 2 ) =1 ll/llAfl(M 2 ) =1 II/IIa/1(B 2 ) =1 

< _ su p II/IIm~(r)II/IIm i (r 2 )||sIIm i (r 2 ) = ll/IU/~(R 2 )llffllM i (R 2 
II/IIm1(e 2 ) =1 



A-/!(R 2 ) 



Note that with <j)*(t) = </)(—t), we have 

r{t)<j>(x - t)e- 2mt ^dt = V^r(x, £), 

which is a bandlimited function since 

^ V r r(x,Cie 27Tlt ^ xl ' dxdS, = J r{t)(t){x ~ t)e- 2mxv dx = r(t)(p(v) e 

Using that the Ai^-norm is invariant under partial Fourier transforms and the equivalence between 
the M 1 and L 1 norms which is implied by the bandlimitation of P/ 1+ [_7v( e )jv( e )](a; + q) V^*r(a;,£) 
to (-l/T,l/T)x{0} + (-6,ft + S)x(-6,T + S), we obtain 



2iriti> 



Ph + l-N(e),N(t)](x + q) r(t) <f>(x - t) 



M!(R 2 ) 



Ph+[-N(e),N(e)](x + q) Vpr(x,£) 
Pl 1+ [-N(e),N(e)]{x + (?) V<j,*r(x,£) 



M 1 (R 2 ) 
)' 
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Fix g G <S(R) and observe that || ^/|| ip(K 2 ) defines a norm on M P (R) equivalent to the M P (R) 
norm given in (|4.1I) For any Acl we obtain the uniform bound 

|| Y C n 6 nT \\ M °o {R) >i\\V g ^ c n S nT\\L°°(R) = \\ Cn d{ nT ~ t)e 27TtmT \\ L c (R) 

nTeA nTeA nTeA 

< II X! I c ™l \9( n T — £)|||l°°(r) < || |c n | l5( n ^ ~ ^)|IU°°(R) < °°- 

The first norm inequality stems from the fact that for all g £ M 1 (R), || V^/H^p^) defines a norm 
on Af p (R) equivalent to the M P (R) norm given in (l4~Tj) . 

Combining this upper bound on ||k; — W2\\m°°(r) with the above estimate for \\<jh c ~ (j h\\l° c {S) 
and (|5.2[) , we conclude 

^ £ , ii-ffciU(M~(R))T^7- hmu°° 

< De||i?c|U(M~(R)) - De\\aH G \\L°°(R2) < De\\aH\\L°°(m?) < De^i. 

Choosing R(e) above large to yield D small enough to compensate all the multiplicative constants, 
we obtain 

As a meaningful statement is only obtained for e < 1, this bound directly implies that 

ll^if II i°°(K 2 ) < 2 M- 

Combining this with the bound 

\Wh - <^h\\l^(R 2 ) < Ikff - Cffclli""^ 2 ) + \\cfh c - o#-||l~(r 2 ) < 2e/i, 
Theorem 12.71 directly yields the result with a constant of twice the size as in Theorem 12.71 □ 

6. Reconstruction of bandlimited operators from discrete measurements 

This section concerns the discrete representation given in Theorem 12.91 First, we prove this 
theorem, hence establishing that indeed this representation is globally exact. 

Proof of Theorem \2.9[ The proof is similar to the proof of Theorem l2. 41 given in (23]. The main 
idea is to use a Jordan domain argument to cover a fixed compact set M of size less than one by 
shifts of a rectangle that still have combined area less than one and then to combine identifiability 
results for each of them to obtain identifiability for the whole set. Indeed, there exist L prime and 
T, O > with TO = \ such that 

L-1 

supp(r?) C (J R + (kjT, njCi) C [—{L - l)T/2, (L + l)T/2] x [-LO/2, LSI/ 2] 
3=0 

= [-1/(20) + T/2, 1/(20) + T/2] x [-1/(2T), 1/(2T)] 

where i? = [0, T) x [—0/2, 0/2), and the sequence (kj,rij) € Z 2 consists of distinct pairs. For S > 
small enough (and possibly slightly smaller T, O, and a larger prime T), one can even achieve 

L-1 

M s C (J R + (k j T,n j Sl) C [-(T-l)T/2,(L + l)T/2] x [-LO/2,LO/2] 
j=o 

where Ma is the ^-neighborhood of M. 

Fix such 5 and let r, € 5(E) satisfy (|2~H)) and (|2~15)) for this <5. Clearly, 

(6.1) (k,n) ^ (k 3 , nj ) for all j implies S 5 n (i?+(JfeT, nSlfj = and r](t,-y)r[t-kT)^-nS}) = 0, 

a fact that we shall use below. 
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h + l 



-N(e),N(e)](x + q)V r r(x,£) 



TM 



Define the identifier w = X) n ez c « ^nT> where {c„} is L-periodic and observe that 

= If v(t,l)e 2 "' /{x - t} Y, c kSkT(x-t)dtd 1 
= Y. Ck I ri(x-kT, 7 )e 2 ^ kT d7 

L-l . 

r-H 1, — n " 



m€Z fc=0 



for any p e Z. We shall use the non-normalized Zak transform Z^t : L 2 (R) — >• L 2 ([0,LT) x 
[-0/2,0/2)) defined by 

We compute using the Poisson summation formula and the fact that O = 1/LT 
(Z LT o H)w(t, v) 

= Y Hw(t - nLT) e 2ninLTv 

ra£Z 

= J- e ^TnLu J2 Ck+p J rj(t - (nL + mL + k+ p)T, 7 ) e 2 «^ mL + fe +rt T d 7 
= E c *+p E e 2mT " L " /" f|(t -K + H p)T, 7 ) e 2 ^ T ((™-») i+fc +rtd7 

fc— m.ntEZ 
L-l „ 

= E Cfe +f E v(t~(mL + k+ p)T, 7) e 2 -^( mi + fc +rt T £ e 2-ni(,-7)T d7 

fc=0 meZ nSZ 

= E c ^+p E [ r](t~(mL + k+ p)T, 7) e 2 -^+ fc +^-±. £ _ 7 )rf 7 

/c=o mez nez 

L-l 

= ^ C fc+p Y n(t - {mL + k+ p)T, V + nfl) e ^i{^n){mL+k+p)T 

By (j6TT]) we get for p = 0, . . . , L - 1, 
(6.2) r(t)^(v){Z LT oH)w{t+pT,v) 

L-l 

= 0^c fc+p ^ r(t)$(v)ri(t - {mL + k)T,v + n9)e 2 * lT ^ +nn ^ mL+k+ rt 

k—0 m,nGZ 
L-l 

= Q ^ cp +kj r(t)<j»(v)ri(t + kjT, v + n,-fi) ^.(f+n^ro+k^ 

.7=0 

L-l 

= n e 2nwpT ^(r^M"^c) p (e 27rw ^ T r(i)^) r;(i + fyT, 1/ + njfi)) , 
i=o 

where here and in the following, T : (cq, c\, . . . , cl-2i cl-i) H> (cl-i, cq, . . . , cls, Ci_a) and 
M : (c ,ci,...,c L _ 2 ,c L _i) k> (e 2 " / L co, e 27ril / L ci,..., e 27ri ( L " 2 )/ L c L _2, e 27ri ( i ~ 1 )/ L c L _i), that 
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is, (T kj M nj c) p — e 27Tl 3 l 3 Cp+kj- Equivalently, we obtain the matrix equation 

(6.3) [e~ 2 ^ T r{t)^{v){Z LT o H)w(t+pT, v )]f^ 

where A is a L x L matrix, whose jth column is T kj M Uj c E C L . A is a submatrix of the L x L 2 
marix G, whose columns are {T k M l c}^jJ^ . It was shown in [TO, that if L is prime, then we can 
choose c € C L such that every L x L submatrix of G is invertible. In fact, the set of such c G C L 
is a dense open subset of C L [16]. Hence we can apply the matrix A -1 —: [bj P ]j p=1 on both sides 
of Equation (|6.3j) to obtain 

(6.4) e 2wwk J T r(t)4>(v)r](t + k^T, v + n^) 

L-l 

= LTY^ b w e- 2 ™ pT r(t)^{v){Z LT o H)w(t + P T, v) 

p=0 

for every j = 0, 1, . . . , L — 1. 

In fact, until this point the proof agrees with the proof of (|2.8I) in Theorem 12.41 Indeed, if we 
extend {bj p } p to a L-periodic sequence by setting &j iP + m z, = bj P , replace the so far unused property 
(|2.15[) by (|2.9j) then further computations [22] give 

L-l 

h(x, t)=LT^2 r(t - k T) ( ^ b jq Hw(t - {k 3 + q)T)4>{x - t + (k 3 + g)T)) tM*"M*-*). 

3=0 qGZ 

Observe that f|2.15[) implies that (r, TZ x ^Z) = {7fcT-A^^Ln/^ 2 r }fc^ez is a tight Gabor frame 
whenever /3 2 > 1+25/Tas, in this case, (r, ||ZxiZ) = (r, ^TZxflLZ) is an orthogonal sequence 
and the Ron-Shen criterion applies [Tl] [26] . The same arguments imply that (<^>,f2Z x -^Z) is a 
tight Gabor frame. Using a simple tensor argument, we obtain that m,n,i,k\m,n,i,k& forms a 
tight Gabor frame where 

*m,rM,fc(*> V ) =T( k T,nQ.)M L (m(3 2 T,T/h) r<S)(f>(t,v) 

= e M ~i + ~2 > r (t - kT) 4>{v - nVt) . 

The frame bound is im 2 ™/(/3i/3 2 ) = VG^iAO- We set <S> m -n,l-k = Fs^m^k- Clearly, as 
!F S is unitary, we have that {^> m ,n,i,k}m,n,i,k<£i forms a tight frame with frame bound l/(/3i/?2), in 
fact, a tight Gabor frame as 

= e 2 ^ nm+kl ^ x (T lL n/p 2 M kT m) (T m TL /f3l M n ^)(x). 

Note that (|6.1[) together with the fact that the symplectic Fourier transform is unitary implies 
that the coefficients in the Gabor frame expansion of <x satisfy 

(a^m-n^i-kj) = {ri,^m,n,l,k) = unless (n,k) = (rij,kj) for some j. 
Hence we need to estimate a~^ i = (a, ^ m ^ njt i-k ,-) f° r j = 0, 1, . . . , L — 1. We obtain by (|6.4[) 
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a m,t = ( (T ^m,-n J ,L-k j ) = (r), *m,nj ,l,kj ) 



= / / 17(4, ^)e~ 27 " L( 31 + — K — > r (t - kjT) <$>{y - n^dtdu 

= JJ r(t) ^>(u)rj(t + kjT, v + nj^l)e 2vivkiT e7 2 ^ L ^ + ^ )+v ^ T dtdv 

= [[ LTJ2 b Jp e- 2mvpT r{t)^{ V ){Z LT o H)w(t + pT, v) ^(H^+^+^t) ^ 

= LT J2 bj P [ [ r{t)${v) e- 2 ^ T {Z LT o H)w(t + pT, v) ^iH^+^+^T) ^ 
= LTj2b jp ff r(t)^{v) e- 2 ^ T J2Hw(t+ P T-qLT) e 2 ™« LT e - 2m(L( ^ + ^ )+ ^ T dtdv 
= LTJ2 b j P zZ ( / + P T - qLT)e~ 2niL ^ dt) ( / e 2^-T( g L-p-fe,-mi,/ / 3 1 ) 

= LTj2B jq ( f r{t)Hw{t + qT)e- 2mLi % dt)H e 2 ™ n- q - kj -mL/ 01 ) ^ 

= LTY^BjqtHTi-q-kj-mL/p^f f Hwifye' 2 ™^^ r{t - qT)dt} 
= LT^B jq 4>(T(-q - k 3 - mL/fa)) (Hw,T qT M iLn/h r), 



where Bj g = bj q > for g = rah + q' with g' = 0, 1, . . . , L — 1. 



Set 



Cq.i(Hw) = (Hw,Tq T M eLn/fj2 r). 



In sum, 



1 

(6.5) = L f e -^n +ikjT) £ ,C0 ^ _ ^) ^ _ 0™) 



where 



a%] e = J2 B jq <f>(a(-q - ki - mL/fa)) C q>l {Hw). 
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Applying the symplectic Fourier transform to (I6.5[) yields 
r,{t,v)=e- m r,{t,v) 



TT 



1 XP2 j=0 m ,££Z 

= e ^ Vt 1^ E E ff S ^T,-„,n) M^.^ r®£(t,i/) 

1 11 2 j=0 m ,£GZ 

LT \ - \ > i, ',>-;.. (,/• r . , / ^ Jtt/'/'/ 



For £/(t, i/) = r®(j> (t, v) e - 27Tivt : we have 

J r s U(x,£) = lj r{t)^{v)e- 2mvt e- 2 ^ t - l ' x) dvdt 

r{t)<j){x - i)e~ 2Tr '^* dt= [ r{t)4>{t - x)e' 27Tl ^ dt = V^r(x, £), 



where we used that real valued implies <j)(y) = <f>{—y). Now, we compute 
a(x,£) =J r s rj(x,£) 



LT 11 



pip2 j=0 m,£eZ 

LT E E ^e 2 ™^ T 7W ( _„.n,_ fcj T)^^_ fe . ri ^_ n . n) ^r(i,0, 



(6.6) = i^y e -2«(xn j a + ^. r)e 2^n fcj T y VflSr(a ._!!^ + A . Tj + 

The convergence in ()6.5[) and ()6.6)) is defined in the weak sense, but can be shown to converge 
absolutely and uniformly on compact subsets. □ 

Next we prove Theorem l2.10l that is, the direct local correspondence between the discretization 
values and the operator action. 

Proof of Theorem \2.1(A We intend to apply Theorems 12.71 and 12.91 We assume that the set 
M as well as its enclosing rectangular grid are fixed, hence also the parameters T, ft, and L. 
The dependence of the constants, auxiliary functions, etc., in the following derivations on these 
parameters will be suppressed for notational convenience; this should be seen as analogue to the 
one-dimensional scenario where the arising constants also depend on the shape and not just the 
size of the frequency support. Furthermore, set Q = max(LT, Lf2). 

We can bound using (|3.3[) 

Wm,e\ = Ko-,$ m _ njA _ fej )| 

(67) < || f || oo 11$™,-™^,-% ||l 

< B|MU|f ®0||i 

For the second inequality, we used that the Li-norm is invariant under translations and modula- 
tions. 
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Furthermore, note that V#r E S(M. 2 ), so there is a decreasing positive function p G <S([0, oo)) 
such that for p{x,^) — p(\x\\ P\\€\) one nas \V<f> r \ — pointwise. 
Now observe that, as p is decreasing, 



CL'l 

OO OO jl 



We use this estimate to bound for arbitrary (x, £) 

kM)l 



L-l 



PlP2 J=0 ( m LT//3i,lLn/p 2 )eS Pl P2 



<- 



<- 



E i 

(mLT/Px,lLn/P2)eS 
L-l 

E E 



T Ci) i 



Pl P2 



: A/32 ^— ' ^ 8CT 



( ,mL 



k 3 )T 



,£L \ s 

e _(_ + ni )n) 



L — 1 oo 



M , /mL\ LQ (IL 

<f (IIpIIi + IIpIIoo) 



=: C lM . 



and hence 

Ik -cf||oo< IkHoo + ll^lloc <M + f(lH|l 

By the definition of 5, for every S > 0, there is a constant C(<$) such that for any fixed < j < L, 

(6-8) 5 >8||p|| L i( R+) ||p|| L i[ C ( 5 )_2Q,oo) > 8\\p\\L 1 (([-C(5)+Q,C(S)-QYY), 

and hence, for (x,£) e S — B(C(S)), 

L 1 y-\ / / //</, 



(6.9) 



5 >- 



/3l/?2 



V /3i 



■kj)T 



/IL 



To obtain (|6.9p from (|6.8p , the boundary term in the discretization of the integral and the shifts 
by kj and rij, respectively, are each compensated by increasing the dimensions of the integra- 
tion/summation domain by LT and Ltt in time and frequency, respectively, both of which are 
bounded by Q. 

Note furthermore that, as (x,£) E S — B(C(5)), a necessary condition for 



-C(S),C(S)] 



is that 

(mLT/PiJLQ/fo) i S. 
Thus, using (I6.7[) and the triangle inequality, we can bound (|6.9p from below obtaining 



(6.10) 



S/j, > 



L 2 T 



E 



(mLT/l3 1 ,lLQ./p2)<tS 



9. 



2:S 



Hence forming a weighted average (with complex weighting factors of modulus one) of Equa- 
tion (|6.10p over the L choices of j , we obtain 



5n > 



LT L ^ 2m(xnj n + t kj T )e 2. mj n kj T £ o^r(x - (^ + **)t, £ - (~ 



-n 



in 



This yields \\a — ^\\l x (S-b(C{5))) < ^M- Hence by Theorem 12.71 we conclude that 

\\Hf-Hf\\2<C^njL 



for all functions / which are ^4— time-frequency-localized to S — B(C{5)) — B{d{e)). The result 
follows by choosing 8 = min ^^=r£, C\ej and Z?(e) = C(<5) + d(e). □ 

7. Quantization of bandlimited operators 

The underlying idea of the quantization schemes in Theorem 12 . 1 1 1 and Corollary |2.12l is based on 
EA modulation. At the core of this arguably most influential coarse quantization paradigm is the 
observation that the reconstruction formula in Theorem 12 . 1 1 directly corresponds to the application 
of a low-pass filter. The key idea is then that, while, due to the coarseness of the alphabet, the 
sequence y n — q n cannot be made uniformly small, it can be chosen to be approximately high-pass, 
that is, close to the kernel of the low-pass operator; hence it almost vanishes in the reconstruction 
procedure. In other words, f(t) = T^2 n£Z q n (j){t — nT) is a good approximation for /. 

The main goal of this section is to show, in general terms, the possibility of combining the 
discretization procedure presented above with coarse quantization schemes. Hence, we will restrict 
ourselves to the simplest possible EA modulator, a so-called first order SA modulator. Usually, 
with so-called higher order modulators [51 [7] considerably better error decay rates can be 
achieved. There are no specific obstacles that would prevent the direct application of such higher 
order modulators in the operator context. However, the estimates would be considerably more 
complicated without providing much additional insight, which is why we we refrain from presenting 
them here. 

To define a first-order one-bit EA modulator, we first fix a kernel (j) such that </> € S and 



supp (j) C 



An Ap 
2 ' 2 



for some Ao > fi. Then the minimal sampling step is 7o = j^- So sampling 

at step size T corresponds to an oversampling ratio of A = As we are interested in the 
approximation behavior when the redundancy of the dictionary increases, we do not want to 
consider multiple quantization schemes. Hence we fix an underlying quantization alphabet with 
an associated quantization rule. Again, we will focus on the simplest possible scenario, namely 
one-bit quantization, where the quantization alphabet just has two elements. 

Using this quantization rule, a first order EA modulator computes a sequence of quantized 
values {q n } by means of the iterative scheme 

(7.1) u n = u n -i + f{nT) - q n 

q n = sign(/ (nT) + 

with an initial condition uq. 

One can show (see for example [3]) that if \f(t)\ < 1 and |uo| < 1, then the state variable u in 
(|7.ip satisfies 

Mn e 1, \u n \ < 1. 
This entails a bound for the quantization error of 

1 



11/ - /lli^(R) < 
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Proof of Theorem \2.11\ By Theorem l2.91 a = LT Y]j—i uijcr^, where uij are phase factors and 

PlP2 ^„ Pi P2 

We will quantize the doubly- indexed coefficient sequence {cr^A of the expansion in (|7.2p , sepa- 
rately for each j = 0,l,...,L — 1. For simplicity, we will drop the superscript (j) for now, since 
most of the proof proceeds independently for each j. We also drop the rij and kj, as for fixed j, 
they can be absorbed in x and £. Hence we seek to quantize an expansion of the form 

We proceed by applying first order Sigma-Delta modulators subsequently in both time and fre- 
quency. Such an approach has been successfully applied by Yilmaz [2 7) to devise a coarse quanti- 
zation scheme for Gabor expansions. 

More specifically, we recursively define doubly indexed sequences u,p, v, r as follows. 

Pni,ri2 — sign(ti ni _i ;7l2 -t-fj ni n2 ). 

(7.2) 



sign(u ni , n2 _i + u ni ,„ 2 ). 



For a bivariate sequence a, we denote (Aio)„ 1)Jl2 = a ni:1l2 — a ni _i iT[2 and (A 2 a) ni n2 = a ni n2 — 
a ni.n 2 -i- Then we have 

(AiAauJm.na = o" n n — g ni „ 2 , 
(7-3) _ . 

Qni,n2 — Pni,ri2 ' v^l^yni ,Ti2 ■ 

The resulting sequence q m ,t G {±1,±3} of quantized values will be used for reconstructing cr, 
these equation hence define the bivariate Sigma-Delta modulator. Note that to properly define the 
Sigma-Delta modulator, one also needs to set an initial condition for some finite index. However, 
the above recurrence relations are reversible in the sense that they allow to apply the quantization 
procedure backwards in time and frequency, hence allowing for the acquisition of quantized values 
for all (m,£) G Z 2 . A more realistic scenario with only finitely many quantized values used 
for the approximation is obtained by combining this result with the localisation result given in 
Theorem ED] 

As the Sigma-Delta scheme given by (|7.2[) and (|7.3I) is a combination of two first order Sigma- 
Delta schemes, its stability is a direct consequence of the stability of (|7.1[) . see also [37] • We 
summarize this observation in the following proposition. 

Proposition 7.1. Suppose that |o" niin2 | < 1. If |uo,n 2 | < 1 f or each 712 G Z, and \v nu o\ < 1 for 
each ri\ G Z, then |um,n 2 | < 1, 1 1*711,712 1 < 1 f or every 711,712 G Z. 

Stability is a crucial ingredient for the error analysis of the bivariate Sigma-Delta scheme. 

Proposition 7.2. Let a : R 2 ->■ C be as in (|7.2[) and g TOi ^ to, I G Z fee f/ie quantized values 
resulting from the bivariate Sigma-Delta quantization given by (|7.2I) and (|7.3p loii/i er m ,£ < 1 /or 
aH to, ^ G Z 2 . TTien i/ie symbol reconstructed according to the formula 

1 / to£„ . £L. 

P2 



satisfies 



1 x / mL k,u 



1 a 2 

0- - cr £,«. < — — - ——V^r . 

P1P2 oycte l 1 



2:, 



Proof. We compute 



= | ^ (A 1 A 2 v) m ,e 

= I E * 



(m+l)LT/Pi q £ 



(m+l)LT/0i f (£+l)LQ/p 



±V>r(x-u,t 



Lit) du 



< 



7 y " m ,{ . . 

m ^ eZ JmLT/Pi J£Ln/f3 2 ? 

r (m+l)LT//3i r (t+l)LSl/fa q2 



V] \ v m,e\ / / 

' 4 ^ JmLT/p 1 J I. 



m,ieZ 



< 



if 



-V^r(u, v) 



d£ dx 



V^r^x — u, £ — v) dvdu 
V 4> r{x -u,£-v) 



dvdu 



dvdu . 



I dy dx 

In the last step we used the stability of the scheme, as established by Proposition 17. II 



□ 



To complete the proof of Theorem l2.11l we recall from (|6.7j) that there exists a constant C such 
that, for all m,£ <G Z, |er^J < CHctHoo < cC by assumption of the Theorem. Hence choosing 
c = 7!) we ensure that |<7„J < 1 and we can apply Proposition 1 7 . 2 1 to conclude that 



<9£ 9a; 



The right hand side is bounded by an absolute constant, so by Proposition l4.1l for p = 2, we conclude 
that \\HWf - H u) f\\ L 2 < a^H/IU^, where and ffO') are the operators with symbol 
and er^ , respectively. Choosing C — L 2 TC, the theorem follows via the triangle inequality. □ 



Sketch of proof of Corollary \2.12\ We quantize H as described above. Due to the recursive 
nature, we need only finitely many of the o~ m ,£ to compute any given finite set of q m ,e's. This 
gives rise to the first summand in the error bound. The step from the quantized representation to 
the quantized representation with only finitely summands proceeds completely analogously to the 
proof of Theorem 12. 101 in particular, the error arising in this step can be bounded by the second 
summand of the error bound. Due to the completely analogous nature of the proof, we will not 
repeat the details here. 
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